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ABSTRACT 

Context. Flickering is a ubiquitous phenomenon in cataclysmic variables (CVs). Although the underlying light source is one of the 
main contributors to the optical radiation, the mechanism leading to flickering is not understood as yet. 

Aims. The present study aims to contribute to the set of boundary conditions, defined by observations, which must be met by physical 
models that describe the flickering. In particular, time lags in the occurrence of flickering events at different wavelengths over the 
optical range are examined. 

Methods. To this end, the cross-correlation functions (CCFs) of numerous light curves of a sample of CVs are analysed that were 
observed simultaneously or quasi-simultaneously in different bands of various photometric systems. 

Results. Deviations of the maxima of the CCFs from zero time-shift indicate a dependence of the flickering activity on the wavelength 
in the sense that flickering flares reach their maxima slightly earlier in the blue range than in the red. While the available observa¬ 
tional material does not permit detecting this individually in all observed systems, the ensemble of all data clearly shows this effect. 
Particularly instructive are the cases of V603 Aql and TT Ari, where time lags of 15 s . 1 and 4 S .3, respectively, are observed between the 
U and R bands. In principle this can be understood if during the development of a flickering flare the radiation characteristics of the 
light source responsible for flickering change such that in the early phases of a flare more short-wavelength radiation is emitted, and 
later on, the peak of the emission shifts to the red. Respective scenarios are discussed and shown to be in qualitative and quantitative 
agreement with observations. 
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1. Introduction 

Cataclysmic variables (CVs) are well known to be short-period 
interacting binary systems where a Roche-lobe filling star, the 
secondary, transfers matter via an accretion disk to a white dwarf 
primary. One of the most striking photometric characteristics 
common to all CVs are short-term variations of the optical light 
occurring on time scales of seconds to a few dozen minutes with 
amplitudes of some hundredths of a magnitude up to more than 
an entire magnitude. This phenomenon is called flickering and 
appears to the eye as a continuous sequence of overlapping flares 
in the light curves of virtually all CVs (but note that this does 
not necessarily mean that flickering consists of a linear superpo¬ 
sition of independent events). Indeed, it can safely be said that if 
it does not flicker, is is no CV. 

Flickering has been observed for more then one and a half 
centuries: In 1856 Pogson (i 18571 1 visually observed variations 
of the dwarf nova U Gem that might be interpreted as flicker¬ 
ing. So did Baxendell in 1858 (Turner, [1907 ). Modern obser¬ 
vations of flickering started in the middle of the twentieth cen¬ 
tury (Henize 119491 Linnell 119501 Lenouvel & Daguillon 119541 
Johnson et al. [l954l Walker & Herbig [l954l Walker [19571 ). The 
amplitude of the flickering immediately shows that the underly¬ 
ing light source is a major contributor to the optical light of the 
respective systems. On the mean over all CVs, it emits at least 
25% of the radiation of the primary but can easily contribute as 
much as 50% and in strong flares can outshine the non-flickering 
light sources even by a factor of almost 2 (Bruch [T989111992]) . 
Even though flickering is thus due to one of the optically most 


luminous light sources, it remains among the least understood 
phenomena observed in CVs. 

Flickering is not restricted to CVs or to optical wavelengths. 
The phenomenon appears rather to be intimately related to ac¬ 
cretion of matter from a disk onto a central body. Herbst & 
Shevchenko (1 1999k for instance, interpreted flickering-like vari¬ 
ations in pre-main-sequence variables as being due to unsteady 
accretion (see also Kenyon et al. 12000 ). A comprehensive re¬ 
view of aperiodic X-ray variability in X-ray binaries with neu¬ 
tron stars and black holes was presented by van der Klis (12004k 
Flickering-like variability also occurs in active galactic nuclei 
[see e.g. Garcia et al. d 19991) and references therein]. The time 
scale of the flickering-type variations depends on the dimensions 
of the underlying system, and the frequency range where it is 
best observed depends on the compactness of the accretor. 

The first encompassing systematic study of the flickering 
properties in cataclysmic variables has been performed more 
than two decades ago by Bruch (1 989 . 1992k who defined sev¬ 
eral statistical parameters to describe the phenomenon. In addi¬ 
tion to quantifying the already mentioned strong contribution of 
the flickering light source to the total optical light, he showed 
that the increase of its spectrum to shorter wavelength is steeper 
than can be explained by simple models such as a black body 
or a power law. Investigating various physical scenarios, Bruch 
( 19891 [19921 concluded that instabilities in the inner accretion 
disk and/or the boundary layer between the disk and the white 
dwarf are most likely to give rise to the observed phenomena. 

Eater work on flickering includes attempts to identify the ex¬ 
act location where it occurs in CV systems. Already in 1985, 
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Home & Stiening (119851 1 had analysed the variations of RW Tri 
through eclipse and concluded that the flickering light source 
has the same centre as the accretion disk, but is less extended. 
This was later confirmed by Bennie et al. (Il996l >. Similarly, 
Welsh et al. (11996b found that flickering in HT Cas is concen¬ 
trated to the centre of the disk, but is not confined to the im¬ 
mediate vicinity of the white dwarf. Modifying the technique 
introduced by Horne & Stiening ( 119851 ) to avoid systematic ef¬ 
fects, Bruch d 19961 2000 1 found similar results for several more 
eclipsing CVs, but that to a lesser degree the hot spot also con¬ 
tributes. However, in IP Peg the hot-spot flickering dominates. In 
V2051 Oph, Baptista & Bortoletto (120041) identified two flicker¬ 
ing components, namely a slowly varying one that they associ¬ 
ated with the disk overflowing gas stream from the secondary 
star, and a more rapidly varying one with a radial distribution 
equal to the steady disk light. 

The behaviour of flickering in the frequency domain has 
been the subject of many studies. Fritz & Bruch ( 119981 ) stud¬ 
ied wavelet transforms of numerous light curves of many CVs. 
They parametrized the resulting scalegrams (Scargle et al.i l 993) 
in terms of its inclination a (frequency behaviour) and its value 
X (flickering strength) at a reference time scale. For a given sys¬ 
tem, a and X are stable over many years. On average, flickering 
is somewhat bluer on short time-scales than on longer ones. CVs 
of different types and photometric states occupy distinct regions 
in the a - E - plane. A similar study limited to the intermedi¬ 
ate polar V709 Cas was published by Tamburini et al. (120091 ). A 
wavelet study in X-rays was performed by Anzolin et al. ( 120101 ). 
The X-ray data are distributed in a much smaller area of the 
a - L parameter space than the optical data. The authors ex¬ 
plain the similarity of the X-ray flickering in objects of different 
classes together with the predominance of a persistent stochastic 
behaviour in terms of magnetically driven accretion processes 
acting in a considerable fraction of the analysed objects. 

Many studies of the frequency characteristics focus on power 
spectra where flickering light curves in general exhibit a be¬ 
haviour similar to red noise, that is, on a double logarithmic 
scale, a linear decrease of power towards high frequencies (but 
with a flattening towards low frequencies; Bruch fl989l [1992; 
Schimpke ll998| . Dobrotka et al. (120141 ) detected two red-noise 
and two white-noise components in the X-ray power spectrum 
of RU Peg, indicating the presence of two turbulent regions. The 
long and low-noise high-cadence light curves of some CVs ob¬ 
served by the Kepler satellite represent superb data from which 
to study flickering. Scaringi et al. ( 2012al ; 20l2bl) used these 
data to construct the flickering power spectrum of MV Lyr, 
which they decomposed into a broken power law and a series 
of Lorentzian components, following Belloni et al. ( 120021 and 
Novak et al. ( 119991 ). They detected a linear rms - flux relation 
such that the light source becomes more variable when it be¬ 
comes brighter, as has earlier been observed in X-ray binaries 
and AGNs (Uttley & M c Hardv l20m1 Uttley et al. lMBI Heil & 
Vaughan 20101 ). This is interpreted as an indication that flicker¬ 
ing variations are due to a multiplicative as opposed to an addi¬ 
tive process, such as shot noise, since in the latter case the rms 
- flux relation would be destroyed. In addition to MV Lyr, rapid 
variations were also analysed in Kepler data of other CVs by 
Scaringi et al. (120141) . 

Interpretations of the frequency behaviour in terms of the 
propagation of disturbances in the accretion disk leading to 
flickering have repeatedly been discussed in the literature. For 
instance, Lyubarskii (119971 ). Yonehara et al. (119971 ). Pavlidou 
et al. (120011 ). and Dobrotka et al. (120101 ) developed conceptu¬ 
ally similar models that differ in details concerning the origin 


and the propagation of the disturbances. Of these, the model of 
Lyubarskii (119971 ). formulated originally to explain X-ray varia¬ 
tions in galactic and extragalactic X-ray source, and which in¬ 
vokes viscosity fluctuations at different radii, has found most 
widespread attention. In the optical range it has sucessfully been 
applied to MV Lyr by Scaringi (120141 ) [but note that Dobrotka 
et al. (120151 ) was able to interpret the same observational data 
also in the framework of the model of Dobrotka et al. (120101 )1. 
This scanario is quite attractive, in particular because it seeks to 
explain optical flickering in CVs and X-ray fluctuations in neu¬ 
tron star and black hole binaries as well as AGNs in a unified 
model. However, it has not yet been convincingly shown that it 
can explain all flickering characteristics observed in CVs (or, al¬ 
ternatively, that certain properties cited in the literature are based 
on incorrect assumptions or faulty techniques). Therefore, alter¬ 
native models should not yet be discarded. 

None of these mentioned models is concerned with a de¬ 
tailed description of the emission mechanism that leads to the 
observed flickering. Of the few attempts to construct a realis¬ 
tic physical model for the temporal and spectral evolution of a 
flickering flare, the work of Pearson et al. (l2005t) deserves to be 
mentioned, which I discuss in more detail in Sect. 15.21 

As a broad-band phenomenon, in the optical range flicker¬ 
ing is apparent from the blue to the red atmospheric cutoff and, 
not surprisingly, in different photometric bands it is not inde¬ 
pendent from each other but strongly correlated (see, e.g. Fig. 2 
of Bruch 1992). However, depending on the temporal evolution 
of the properties of the emitting light source (e.g. temperature 
and opacity), it is conceivable that individual flares do not reach 
their highest intensity simultaneously at all wavelengths, thus 
leading to time lags between the flickering in different photomet¬ 
ric bands. Detection and quantification of such lags will provide 
constraints on physical mechanisms explaining flickering. 

A temporal shift of the flickering at different wavelengths 
can be traced by measuring the location of the peak of the cross¬ 
correlation function (CCF) of light curves in various bands. This 
technique has been applied by Jensen et al. ( 119831 ). who observed 
in this way a lag of about l m between the flickering in the opti¬ 
cal and X-ray bands in TT Ari. Similarly, Balman & Revnivtsev 
(120121 ) detected lags between 96 s and 181 s between UV and X- 
ray observations in five more CVs. Using optical observations of 
AM Her, Szkody & Margon (119801 ) deduced from the asymmetry 
of the maximum of the CCF a lag of the flickering between the 
U and B bands of 8 s in the sense that variations appear earlier 
in U than in B. On the other hand, Bachev et al. (120111 ) found 
no conclusive evidence for delays of the flickering in different 
bands of the nova-like variable KR Aur. This may be due to the 
comparatively coarse time resolution of 30 s of their data. Bruch 
( 1989 , 1992 *) also obtained null results, but this can probably be 
attributed to insufficient quality and quantity of the observational 
data. 

A more general approach makes use of the coherence func¬ 
tion, that is, a Fourier frequency-dependent measure of the linear 
correlation between time series measured simultaneously in two 
spectral ranges (Vaughan & Novak IT99 7). This permits detect¬ 
ing time lags as a function of the (temporal) frequency of the 
underlying signals. This technique has found widespread appli¬ 
cation in the study of X-ray binaries and AGNs, but has only re¬ 
cently successfully been applied to optical data of CVs. Scaringi 
et al. (12013j> studied the Fourier-frequency-dependent coherence 
of light curves of MV Lyr and LU Cam observed at high time 
resolution (0 S .844 - 2 S .276) in the SDSS bands u', g' and r'. 
They also described in detail the applied data reduction method. 
Scaringi et al. (120131 ) found that low-frequency variations (be- 
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low =3 10~ 3 Hz) occur about 10 seconds later in the in r' than 
in u' in LU Cam. The corresponding time lag is « 3 seconds in 
MV Lyr. At higher frequencies, the errors become too large to 
permit reliable measurement. 

Here, I investigate numerous light curves of various CVs ob¬ 
served in several photometric systems to measure time shifts of 
the flickering as a function of wavelength. As mentioned, the 
time lags observed by Scaringi et al. (120131 1 are most signifi¬ 
cant for low-frequency variations. Moreover, the (approximate) 
power-law behaviour of the flickering power as a function of 
frequency (Bruch 11989 '. 199 2) means that low-frequency varia¬ 
tions strongly dominate high-frequency flickering. Considering 
furthermore that the CCF will be most sensitive to the strongest 
variations as well as the limitations of the available data (time 
resolution, noise, length of individual data sets; see also Sect. 0, 
I prefer here to use the cross-correlation function, which is prob¬ 
ably more robust than the much more sophisticated method ap¬ 
plied by Scaringi et al. (i2013t >. Not using the coherence function 
implies the loss of some information. But in view of the data 
limitation, this disadvantage is compensated by more reliable 
results than would be possible using more sophisticated tech¬ 
niques. To measure the time lag of the flickering as a function 
of wavelength, I therefore examine the CCFs of the light curves, 
that is, the normalized cross covariant as a function of the shift 
of the data in a reference band with respect to the other (compar¬ 
ison) bands of the photometric system. In Sect. [2] the reduction 
method is explained. The observational data are briefly intro¬ 
duced in Sect. [3] before the results of this study are presented in 
Sect. 0 Scenarios to explain the results are discussed in Sect. 0 
and finally, conclusions are drawn in Sect. 0 


2. Reduction method 


The CCF is calculated in the following way: A light curve is 
expected to have been sampled simultaneously in the reference 
and the comparison band in discrete time intervals separated by 
At. In the data sets used in this study this is often the case, but 
sometimes only approximately so. Then, the light curves are re¬ 
sampled by linear interpolation between data points at the in¬ 
stances t\ + j At, where t\ is the instance of the first data point in 
the reference band, j is an integer number, and At is chosen to 
be the average time interval between data points in the original 
light curve. Comparison and reference data (designated x and v, 
respectively) are then normalized to their average value before 
the mean is subtracted, yielding n pairs of data points (x,y). The 
CCF is evaluated at discrete time shifts r>, such that r/ f = k At 
and k is an integral number. 


ccF(n) = 


koto) 


m^i -1 


v-' m 

x i Z,=i yt+k 


zr= 
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(g=. *<f 


y m 2 _ 1V ,'w) 

1 X i+k m 


Here, i = 1 ...m indicate those data pairs for which reference 
and comparison data overlap after the comparison data have 
been shifted in time by t>. The CCF is thus discretely sampled 
at a time resolution of At. 

If reference and comparison data exhibit similar variations 
but the pattern of variability is shifted by a time interval t s be¬ 
tween them, the CCF will assume a maximum at t = t s = A7 nlax . 

In principal, it is straightforward to determine the maximum 
of the CCF of the data in a reference and a comparison band 
of light curves observed simultaneously or quasi-simultaneously 
at different wavelengths. However, there are a few details aris¬ 
ing from the structure of the individual data sets available for 
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Fig. 1. Example of a light curve and the CCF of the brightness 
variations in two photometric bands. Top: fi-band light curve of 
V603 Aql observed on 1983 July 18. The solid line is a fit of 
a high-order polynomial to the data. Centre: Residuals between 
the original light curve and the polynomial fit. Bottom left: CCF 
of the light curve shown in the top panel with the simultane¬ 
ously observed light curve in the [/-band (broken line) and of 
the residuals between the original data and the polynomial fit 
(solid line). Bottom right: Enlarged section of the maximum of 
the CCF shown as a solid line in the lower left panel (dots) to¬ 
gether with a fit of a high-order polynomial to the data (bro¬ 
ken curve). The broken vertical line indicates the location of the 
maximum of the fitted polynomial. 


this study that deserve to be mentioned. As an example, the up¬ 
per panel of Fig. Q] shows the B -band light curve of the old 
nova V603 Aql observed in the Walraven photometric system 
on 1983, July 18 (see Walraven & Walraven, [19601 and Lub & 
Pel, 1977, for the peculiarities of the Walraven system). In ad¬ 
dition to flickering, it shows another feature frequently observed 
in CVs, that is, variations on longer time scales that may not be 
immediately associated with flickering (at least not in the light 
curve discussed here; see Bruch 1991 ). 

The CCF of the light curves in the B band and the simulta¬ 
neously observed U band is shown as a broken line in the lower 
left frame of Fig. Q] The central peak is caused by the corre¬ 
lated flickering activity in the two bands, but on the whole, the 
profile is rather shallow. This is caused by the longer time-scale 
modulations. Even if these (in a more general case than for the 
present example) were a manifestation of flickering, the limited 
length of the light curves - being of the same order of magnitude 
- would inhibit a reliable measurement of an associated time lag 
of the order of a few seconds between different bands. Therefore, 
to better isolate the variability due to flickering on time scales 
where a lag can be measured more reliably, it is appropriate to 
subtract the slow modulations. This is done by fitting a polyno¬ 
mial of suitable order to the original data in both bands, as shown 
by the solid line in the upper frame of Fig. [I] In all cases when 
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a polynomial was subtracted, the order was chosen such that the 
polynomial followed the long-term variations satisfactorily, and 
the same order was adopted for the reference and the comparison 
bands. The light curve in the central fram of Fig.|T|constitutes the 
residuals between the original data and the fitted polynomial in 
the B band. The solid graph in the lower left frame represents 
the CCF of these residual light curves in B and U. Here, the cen¬ 
tral peak due to the correlated flickering stands out much more 
clearly. 

The lower right frame of Fig.|T|contains a magnified version 
of the upper part of the CCF maximum (solid dots). The coarse 
resolution reflects the rather low time resolution of 21 s of the 
example light curve. To estimate the location of the maximum, 
a polynomial was fit to the data of the CCF Its order depends 
on the number of available points in the peak of the CCF (range: 
-120 s < r < 120 s in this case). This is the dashed line in the fig¬ 
ure. In light curves with a very high time resolution, the CCF of¬ 
ten contains a sharp spike centred on zero time-shift. It is caused 
by correlated atmospheric noise in the reference and comparison 
bands and is masked in the fit process. The quality of the fit indi¬ 
cates that the maximum of the CCF can be determined with high 
precision by calculating the maximum of the fit polynomial. Its 
location is marked by the broken vertical line. The correspond¬ 
ing ATmax is slightly negative in this example. 

To obtain a notion of the accuracy to which A7 m;lx can be 
measured, a simple shot-noise model was used to generate arti¬ 
ficial flickering light curves. While this may not reproduce the 
flickering properties observed at least in some CVs and in many 
X-ray binaries such as the linear rms - flux relation, the depen¬ 
dence of the accuray of A7’ rnax on properties such as statistical 
noise in the data or their time resolution will not depend strongly 
on such details. Thus, 100 artificial light curves were generated, 
each of which subsequently served as reference data. One of 
them is shown in the upper frame of Fig. [2] The light curves 
were constructed in the following way: 1 500 simulated flares 
were distributed randomly on a time base of three hours, which 
is the typical duration of the real light curves studied here. Their 
amplitudes are distributed between 0 m < A < 0™ 1 such that the 
probability of a flare with a particular amplitude to occur decays 
linearly from the minimum to the maximum value of A and the 
probability for a flare with the maximum amplitude goes to zero. 
The flares are symmetric and rise and decay at a mean rate of 
0 m 5 per hour, allowing for a random scatter equally distributed 
between 0 < dm/dt < 1 magnitudes per hour around this value. 
The superposition of these flares, sampled in I s intervals, plus an 
arbitrary constant, constitutes the light curves. To make the light 
curve displayed in Fig. [2] look more natural, random Gaussian 
noise with cr = 0 m 01 was added. Second versions of the same 
data sets were created by applying a constant time shift of 5 s 
to the originals. These served subsequently as comparison light 
curves. Using a constant is acceptable here because a possible 
dependence of the time lag on flickering frequency as observed 
in MV Lyr and LU Cam (Scaringi et al. 1201 3) will not have a 
strong bearing on the accuracy of the A7’ ma x determination, the 
more so because the CCF is - as mentioned - most sensitive to 
the stronger low-frequency variations. 

To study the accuracy of A7’ max on the noise level, Gaussian 
noise was added to the reference and the comparison light curves 
before A7' rnax was calculated. This was repeated for different 
noise amplitudes. The results are shown in the lower left frame 
of Fig. 0 where for each noise level the mean value of A7’ max 
derived from the 100 trial light curves is plotted. The error bars 
are the standard deviation cr around the mean and thus indicate 
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Fig. 2. Top: Artificial flickering light curve used to estimate the 
accuracy of the determination of flickering time lags. Bottom 
left: Dependence of Ar max on the noise level. Bottom right: 
Expected error of A7 , max as a function of the time resolution of 
the light curve for low (0‘”02; open symbols) and moderately 
high (0 m 05; filled symbols) noise levels. 


the accuracy with which A7 , max can be measured in an individual 
light curve. 

The dependence of the accuracy on the time resolution of 
the data was determined in a similar way: The artificial light 
curves were binned to simulate a lower resolution before noise 
was added. Then, the standard deviation of A7’ max was calcu¬ 
lated. It is shown in the bottom right frame of Fig. [2] as a function 
of the time resolution for low (0™02; open circles) and interme¬ 
diate (0 m 05; filled circles) noise levels. 

The accuracy of the results also depends on the length of 
the light curves (or, equivalently, on the number of data points). 
Tests with 100 data sets with half (twice) the length of the origi¬ 
nal ones but with otherwise identical characteristics, subjected to 
the same noise levels as shown in the lower left frame of Fig. [2] 
resulted in errors of Ar max that were 1.44+0.21 (0.70+0.06) 
larger. Thus, not unexpectedly, the errors scale to a high degree 
of accuracy with the inverse of the square root of the length of 
the light curves. 

This exercise shows that AT^x can be measured to an accu¬ 
racy of the order of a tenth of a second in high-quality (long, low 
noise, high time resolution) light curves, but is limited to several 
seconds if they are shorter, noisier and/or of low time resolution. 


3. Data 

A total of 197 light curves of 23 different CVs were used in this 
study. They were observed in several photometric systems over 
an interval of about two decades at different telescopes by many 
observers with a variety of scientific purposes in mind. They are 
thus quite heterogeneous. The data can roughly be grouped and 
characterized as follows: 

VBLUW (Walraven) light curves: They were obtained by 
various observers at the Dutch telescope at ESO/La Silla, 
equipped with the Walraven photometer (Walraven & Walraven, 

11960 ). The instrument permits only a limited choice of configu¬ 
rations, resulting in a time resolution of about 21 s of the present 
data, which may be at the limit of usefulness for the current pur¬ 
pose. The five passbands with decreasing wavelengths from V in 
the visual to W in the ultraviolet are observed simultaneously. 

UBV light curves: Data in this photometric system were gen¬ 
erally observed sequentially in the three bands. This also inhib- 
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Table 1. Mean flickering time lags (in seconds) in cataclysmic variables 


Name 

B/W 

B/U 

B/L 

B/V 

B/R 

B/I 


RX And 




-3.12 

( 1) 




-5.19 

( 1) 





AE Aqr 

-2.16 

( 

1) 

2.59± 

7.97 (55) 

-5.93+17.29 ( 

2) 

-2.82+11.14 (54) 

4.38+ 

7.25 (47) 

5.02+ 

9.77 (31) 

V603 Aql 

-2.40± 

3.53 ( 

4) 

-3.85± 

2.46 (15) 

-3.75+ 

4.18 ( 

4) 

2.44+ 

4.04(15) 

10.58+ 

2.92 ( 7) 

17.15+ 

7.21 ( 2) 

TT Ari 

0.02± 

1.29 ( 

3) 

-1.45± 

1.16(24) 

0.55+ 

0.86 ( 

3) 

1.20+ 

1.30 (24) 

2.25+ 

0.53 (13) 



QU Car 




-1.65± 0.82 ( 3) 




4.90+ 

2.56 ( 3) 





BVCen 

1.88± 

3.04 ( 

5) 

-1.53+ 

1.72 ( 5) 

-3.18+ 

4.12 ( 

5) 

-0.05+ 

4.41 ( 5) 





WWCet 

-1.26± 

2.13 ( 

6) 

1.04± 

3.53 ( 7) 

2.07+ 

4.60 ( 

6) 

0.84+ 

3.47 ( 7) 

-1.78 

( 1) 

-1.77 

( 1) 

TCrB 




0.68± 

1.71 ( 2) 




0.43+ 

3.11 ( 2) 





SS Cyg 




1.06± 

3.12 ( 7) 




4.22+ 

1.41 ( 6) 





HR Del 




-0.75+17.79(19) 




4.23+10.70(19) 

5.39+ 

4.14 ( 4) 



DO Dra 




1.09+ 

0.06 ( 2) 




2.13+ 

0.42 ( 2) 

1.65+ 

0.22 ( 2) 



V795 Her 




1.63+ 

0.01 ( 2) 




0.45+ 

0.58 ( 2) 

-0.07+ 

0.00 ( 2) 



RS Oph 

10.93 

( 

1) 

5.10± 

5.28 ( 8) 

-3.28+ 

0.19 ( 

2) 

3.96+ 

3.29 ( 8) 

8.29+ 

0.97 ( 2) 

10.20+ 

1.37 ( 2) 

V426 Oph 

2.87± 

1.96 ( 

3) 

0.67± 

2.37 ( 3) 

1.63+ 

0.62 ( 

3) 

0.67+ 

1.31 ( 3) 





CN Ori 




-5.64 

( 1) 




0.94 

( 1) 

4.93 

( 1) 



RUPeg 




2.61 

( 1) 




5.28 

( 1) 





IP Peg 




-1.15 

( 1) 




-1.17+ 

0.57 ( 2) 

-0.08+ 

0.49 ( 2) 



GK Per 




0.19± 

2.66(15) 




1.97+ 

2.93(15) 

2.46+ 

2.86 (12) 



AO Psc 




-0.39± 

0.98 (11) 




0.44+ 

1.11(11) 

0.19+ 

0.94(11) 



VY Scl 

-0.15± 

0.83 ( 

4) 

0.65± 

1.21 ( 4) 

0.50+ 

1.90 ( 

4) 

2.66+ 

0.39 ( 4) 





RW Tri 




6.16± 

3.80 ( 4) 




4.35+ 

3.84 ( 4) 

8.81+ 

6.55 ( 4) 



UX UMa 




15.74 

( 1) 




-13.89 

( 1) 

-30.96 

( 1) 



IX Vel 




0.21± 

2.48 ( 5) 




-5.60+ 

8.23 ( 5) 






ited a high time resolution. Most of the UBV light curves re¬ 
garded here are limited to a resolution of ~ 15 s or lower. 

UBVRI light curves: Many suffer from the same disad¬ 
vantages as their UBV counterparts, resulting thus in a rela¬ 
tively coarse time resolution. However, some have been ob¬ 
tained using a photometer equipped with a rapidly rotating filter 
wheel (Jablonski et al. 1994k which permits obtaining quasi- 
simultaneous measurements in all passbands. These data sets 
were observed at a time resolution of 5 s , and are thus more suit¬ 
able for the present purpose. 

UBVR* (Stiening) light curves: These are data sets were ob¬ 
tained with the Stiening photometer (Horne & Stiening 11985 ). 
which is equipped with four filters similar but not identical to 
the standard UBVR filters. All passbands are observed strictly 
simultaneously. The disadvantage that this photometric system 
is not well calibrated is by far outweighed in the present context 
by the high time resolution used to observe the light curves em¬ 
ployed here: Most have a sampling interval of 0 s .5 or I s , while 
some even reach 0 s .2. 

Thus, the most suitable data sets are those observed at high 
time resolution in the UBVRI system and, above all, the Stiening 
data. 

4. Results 

All available light curves were subjected to the reduction method 
outlined in Sect. [2] The B band was used as reference bancQ The 
light curves of all other bands were cross-correlated with the B 
-band light curve, and the time lags AT)^ of the maxima of the 
respective CCFs were determined. In some cases, only parts of 
the light curves were used. This applies in particular to AE Aqr, 
which sometimes exhibits periods of no or very weak flaring ac¬ 
tivity that were then excluded. Eclipses in some light curves of 
eclipsing CVs were also removed. Strong long-term variations 

1 One of the light curves of AE Aqr was observed only in U, V and 
R. In this case, the shortest wavelength band served as reference. 


(time scale: hours), if present, were removed by subtracting a 
fitted polynomial of suitable degree (see Fig. [I]). 

Table [j] contains average values of A7' rnax of the maximum 
of the CCF of each comparison band with respect to the refer¬ 
ence band as indicated in the table header. This is interpreted as 
the time lags of the flickering in these bands. Since short data 
sets should contribute to the average to a lesser degree than long 
ones, the results for the individual light curves were weighted 
with their total duration. The number n of light curves contribut¬ 
ing to the mean are given in brackets. The errors are standard 
deviations derived from all light curves of a given star. No se¬ 
lection was made. All light curves, regardless of the photometric 
system or the time resolution, enter the average. 

4. 1. The ensemble 

No clear picture emerges for the individual systems. In most 
cases, the standard deviation is considerable, often significantly 
larger than the average value of AT m:ix . This is not surprising 
in view of the expected errors determined from the simulations 
performed in Sect.[2]and the limited quality of most light curves. 
Moreover, more often than not (disregarding the errors), there 
is no monotonical dependency of A7 max as a function of wave¬ 
length difference between reference and comparison bands, as 
might be naively expected. However, when the investigated stars 
are considered as an ensemble, at least a trend appears. 

In Fig. [3] (left) the average A7 max for each system is shown 
as a function of the isophotal wavelength /li s0 of the compari¬ 
son band. di so was calculated using the transmission curves as 
listed by Bessell (119901 1 for the UBVRI system and by Lub & 
Pel (1 1 977t ) for the Walraven system [applying a correction based 
on the revision of the mean wavelengths by de Ruiter & Lub 
(119861) 1. together with a first-order approximation of the spec¬ 
tral density S ,\ of a CVs according to the law S ,i oc AT 1 ^, valid 
over a wide wavelength range for an accretion disk radiating lo¬ 
cally as a black body (Lvdnen-Bell ll969 ). The results are listed 
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Table 2. Isophotal wavelengths 


UBVRI system 
Band /lj s0 (A) 

Walraven system 
Band T iso (A) 

U 

3573.2 

W 

3231.9 

B 

4347.6 

U 

3618.2 

V 

5458.6 

L 

3832.6 

R 

6454.9 

B 

4283.1 

I 

7992.0 

V 

5417.9 



3000 4000 5000 6000 7000 80003000 4000 5000 6000 7000 8000 

wavelength (A) wavelength (A) 


Fig. 3. Left: Average A7 max for all investigated CVs as a function 
of isophotal wavelength of the comparison band. Right: Average 
AT max of all individual light curves, weighted in proportion to 
their duration. The broken line is a linear least-squares fit to 
the data points (excluding the points referring to the W and U 
bands), weighted by the inverse of their standard deviations. The 
error bars represent the mean error of the mean. In both frames 
the location of the reference band ( B ) is indicated by a cross. 


in Table [2] Since the transmission curves of the Stiening system 
are unknown, Ti so was taken to be the same as for the respec¬ 
tive bands of the UBVRI system. Moreover, for simplicity in the 
case of the Walraven system the isophotal wavelength of the U, 
B and V bands were substituted by those of the UBVRI system 
in Fig. [3] Finally, A7 max between the B and R band of the only 
light curve of UX UMa was discarded because of a very high 
negative value, deviating strongly from all other results. This is 
due to a peculiar and highly asymmetric CCF. The location of 
the (reference) B band (where Ar max should evidently be zero) 
is indicated by a cross in Fig. [3 

There is a certain trend for Ar max to increase with wave¬ 
length for the photometric bands redward of U, meaning that on 
average there is a slight offset between the flickering in the blue 
and in the visual or red or infrared in the sense that flickering oc¬ 
curs later at longer wavelength. But this trend is not continued to 
the ultraviolet W and U bands. It may not be a coincidence that 
the points deviating from the linear relationship refer to those 
bands that have an isophotal wavelength shorter than the Balmer 
limit (see Sect. IO >. 

The right frame of Fig. [3]is similar to the left one. However, 
instead of showing the average A7' rnax values of the respective 
systems, the average of all individual light curves, weighted in 
proportion to their lengths, is shown as a function of the wave¬ 
length of the comparison band. The error bars represent the mean 
error of the mean (i.e. the cr/ a/h, where cr is the standard devia¬ 
tion and n is the number of data points contributing to the aver¬ 
age). The broken line is a linear least-squares fit (excluding the 
W and U bands) where the data points were weighted according 
to the inverse of the standard deviation. It has an inclination of 


Table 3. Restricted results for individual systems 




V603 Aql 

TT Ari 

AT m ax (B/U) (sec) 

-4.5+2.3 (7) 

-2.0+0.6 (16) 

AT max (B/V) (sec) 

4.1±1.8 (7) 

1.1+0.8 (16) 

A T max (B/R) 

(sec) 

10.6±2.9 (7) 

2.3+0.5 (13) 

A7’ max (B/I) 

(sec) 

17.2+7.2 (2) 

- 

AAT m .JdA 

(sec/A) 

(5.0+0.3) 10- 3 

(1.5+0.1) 10- 3 


(1.9 + 0.3) 10 2 3 * sec/A0. Again, the cross indicates the location of 
the reference band that did not enter the fit. 

The trends observed in Fig. [3] suggests the consistent pres¬ 
ence of flickering time lags in CVs, but may not yet be consid¬ 
ered conclusive evidence. This is the more so because the data 
points shortward of the Balmer jump deviate from the linear re¬ 
lationship. It may be that well-defined time lags in some system 
with high-quality data are “diluted" by large errors from lower 
quality light-curves when regarding the entire ensemble. This is 
suggested by more convincing results on some individual sys¬ 
tems that I discuss subesequently. 

4.2. Individual systems 

While there is thus a trend for a time lag of the flickering at 
different wavelengths when regarding the whole ensemble of in¬ 
vestigated light curves, this is, as mentioned, in general not evi¬ 
dent when regarding individual stars. There are exceptions, how¬ 
ever. Concentrating on high-quality, low-noise light curves with 
a time resolution of <5 S and using as criteria (i) that a reasonable 
number of light curves are available to render the results statisti¬ 
cally reliable, (ii) that the standard deviation of the Ar max values 
of the individual curves is significantly smaller than the differ¬ 
ence of the average A7 max at short and long wavelengths, and 
(iii) that the average A T max is a monotonic function of wave¬ 
length and has a different sign for comparison bands shorter 
and longer than the reference band, two systems were identi¬ 
fied where A7 max has a convincing dependence on wavelength. 
These are V603 Aql and TT Ari. 

For both stars, these results are based almost exclusively 
on the high time resolution UBVR* light curves. Only for 
V603 Aql, additionally, two UBVRI data sets were used. Tabled 
lists the average A7 max values together with their statistical er¬ 
rors. The number of contributing light curves is given in brack¬ 
ets. Figure |4] (upper frames), which is organized in the same way 
as the right frame of Fig. [3] show the same results graphically. 
Clearly, A7 max strictly obeys a linear relationship with wave¬ 
length A, which, in these cases, continues to the U band. The 
derivative of A7 max with respect to A, that is, the inclination of 
a linear fit to the data points weighted by the inverse of their 
standard deviation, is listed in the last line of Table[3] 

If I relax the criteria formulated in the first paragraph of 
this subsection, permitting that in the bands blue-ward of the 
Balmer jump Ar max deviates from a monotonical relationship 
with wavelength, and including also light curves observed with 
lower time resolution (but requiring instead that at least four 
remaining data points are available). Table |TJ shows that the 
relationship between A7 max and A is also quite significant for 
RS Oph (lower left frame of Fig. [4]) for wavelengths longer than 


2 Weighting alternatively the data points according to the number 

of contributing light curves changes this value only slightly to (2.1 ± 

0.4) 10~ 3 sec/A. 
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wavelength (A) wavelength (A) 

3000 4000 5000 6000 7000 80003000 4000 5000 6000 7000 8000 



3000 4000 5000 6000 7000 8000 

wavelength (A) 


Fig. 4. Average Ar max derived from high time resolution light 
curves of V603 Aql, TT Ari and RS Oph as a function of isopho- 
tal wavelength of the comparison band. The broken lines are lin¬ 
ear least-squares fits to the data points (in the case of RS Oph 
without considering the points referring to the W and U bands), 
weighted by the inverse of their standard deviations. The loca¬ 
tion of the reference band (B) is indicated by a cross. The error 
bars represent the mean error of the mean. 


the Balmer limit. For this system, the derivative tlA7' max /tU = 
(3.6 ± 0.5) 1(T 3 sec/A. 

dAr max /d/l differs for the three stars, being more than three 
times higher in V603 Aql than in TT Ari, with the value for 
RS Oph lying in between. This holds true for the average, but 
also for inclinations derived from the individual light curves. 
The magnitude of the flickering time lags thus depends on the 
particular star, but within a given system it appears to be more 
or less stablcQ. However, this statement being based here on only 
three stars still requires confirmation from high-quality data of 
other CVs. 

Scaringi et al. (120131 ) found that the time lags in the stars 
studied by them depend on the frequency of the fluctuations, 
which are larger for variations occuring on longer time scales 
(and even may change sign on short time-scales). The CCF tech¬ 
nique applied here is not able to distinguish between different 
temporal frequencies. Therefore, the observed time lags should 
represent some kind of weighted mean over all time scales. Since 
the CCF is most sensitive to the strongest variations, the mean is 
expected to be biased to the time lags of the dominating fluctu¬ 
ations. Bruch ( 19891 U992| ) has shown that the many CVs these 
occur on time scales of the order of 2 m - 4 m as measured from 
the width of the central peak of the auto-correlation function of 
their light curves. 

The fact that a time lags between the flickering in different 
photometric bands is clearly measured in the two systems with 
the highest quality data available for this study [plus at least 
one additional system (RS Oph) and also the stars observed by 
Scaringi et al. (120131 )1 gives confidence that the trend observed 
in the entire ensemble is not accidental, but a strong indication 

3 In RS Oph the scatter of dAT max /d/i of individual light curves is 
larger than in V603 Aql and TT Ari, which may be due to the inferior 
quality of the available data. 


that a similar time lag is a universal feature of the flickering in 
CVs. 


5. Discussion 

As mentioned in the introduction, reports about a time lag of 
the flickering in different spectral regions have been published 
earlier, but the scope and quantity of these studies remains quite 
limited. In two older papers and in two more recent studies a 
positive detection is claimed. 

The finding of Szkody & Margon (119801 ) was that in AM Her 
the flickering in the V band lags that in the U band by 8 s is not 
based on the measurement of the peak of the CCF, but on its 
asymmetry. They measured the lag at the half-intensity point. 
But it is by no means obvious how an asymmetry in the peak 
of a CCF is related to a time lag between the correlated func¬ 
tions. Indeed, tests with artificial light curves showed that it is 
not easy to introduce asymmetries in the CCFs by shifting the 
flares in one of the bands or by changing their shapcfl However, 
such asymmetries can be caused by longer time-scale variations 
if these are significantly different in the reference and compar¬ 
ison bands. Therefore, it is not clear how reliable the results of 
Szkody & Margon (119801 ) are. They interpreted their findings in 
the framework of the AM Her star model, which does not apply 
to the CVs studied here. 

More relevant in the present connection is the study of Jensen 
et al. (119831 ). who found a correlations between the flickering in 
X-rays and the optical in TT Ari, the same system where I have 
detected a highly significant dependence of A T mdx on the wave¬ 
length in Sect. 14.21 The X-rays lag the optical variations by about 
a minute. Jensen et al. (119831 ) interpreted their results within 
a model where X-rays are emitted by a corona above the ac¬ 
cretion disk. The correlation of variations in the different bands 
can then be understood as due to the efficiency of acoustic and 
magneto-hydrodynamic transport processes from the disk to the 
corona. However, Jensen et al. (119831 ) did not investigate this 
scenario quantitatively. Balman & Revnivtsev (120121 ) obtained 
similar results for additional CVs, but interpreted them differ¬ 
ently: In the framework of a model of propagating fluctuations 
(e.g., Lvubarskii U9971 ). the lag would be due to the travel time 
of matter from the innermost part of a truncated accretion disk 
to the surface of the white dwarf. 

While Szkody & Margon (119801 ) focused on the special case 
of a strongly magnetic CV, and Jensen et al. d 19831) as well as 
Balman & Revnivtsev (12012b compared flickering in widely dif¬ 
fering wavelength ranges, the study of flickering time delays in 
MV Lyr and LU Cam by Scaringi et al. (120131 ) can be compared 
much more directly to the present one. As mentioned, at low 
flickering frequencies (several times 1 O -2 Hz as read from their 
Fig. 3), they measured a delay of ~ 10 sec in LU Cam and ~ 3 
sec in MV Lyr of the r' band variations with respect to those 
in the u' band. With the average SDSS filter wavelengths as in¬ 
formed on the SDSS web page (and neglecting the small differ¬ 
ence between average and isophotal wavelength), this translates 
into a time lag of 3.8 10 4 5 sec/A and 1.1 10 3 sec/A, respec¬ 
tively. These numbers are quite similar to the time lags observed 
here, strengthening the case for such lags to be ubiquitous in 
CVs. 


4 Even the extreme case where the flares are triangular in the refer¬ 

ence band and saw-tooth-shaped in the comparison band (setting the ris¬ 
ing branch of the original triangle to zero) only introduced a very slight 
asymmetry (but changed the location of the maximum significantly). 
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5.1. A heuristic model 

The time lags observed here in the optical bands are much 
smaller than those seen by Jensen et al. (119831) and Balman 
& Revnivtsev (120121 1 and have the opposite wavelength depen¬ 
dence. This is not surprising because the reasons for the lags in 
the optical bands probably are quite different from those for the 
lag between the optical and X-ray range. Here, I find that varia¬ 
tions in the blue occur slightly earlier than in the red. In principle 
this can be understood if during the development of a flickering 
flare the radiation characteristics of the underlying light source 
change such that in the early phases of the flare more short wave¬ 
length radiation is emitted, and later on, the peak of the emission 
shifts to the red. 

To verify that this leads to results that are compatible with 
the observations, a very simple scenario was investigated. I stress 
that this is not meant as a physically realistic model but just to 
show that based on simple but sensible assumptions, it is possi¬ 
ble to reproduce the results found in Sect. [4] 

For this purpose the light source underlying a flickering flare 
was approximated by a black body with a temperature GO that 
over time t varies according to some function 0(f)- The emitting 
area Ch of the black body was also considered a function of time: 
Jt(t). The radiation emitted by the black body as a function of 
wavelength A is then given by the Planck function B\{Q(t)] and 
Jt(t). The radiation flux Fit) detected in a passband of a photo¬ 
metric system is thus 

ph 

F{t) = const B,\ [0(f)] Jt(t) T(A) d/t, 

Ja, 

where 7~(A) is the transmission function of the photometric band, 
and the integration extends between its upper and lower wave¬ 
length cutoff. The constant depends on the distance of the radia¬ 
tion source and is of no relevance in the present context. 

In the absence of a specific physical model for the flickering 
flares, the functions ©(f) and IB(t) are free parameters. They are 
chosen here by trial and error such that the resulting flare pro¬ 
files and their time lag A T msa as a function of wavelength are 
compatible with the observations. Of course, the choices made 
here can by no means be considered as unique. The combination 
of different functions may well lead to similar results. 

The time development of the emitting area is taken to be 
Jt(t) oc f 2 . This is based on the simple notion of a light source 
that expands linearly with time (which may not necessarily be re¬ 
alistic!). With respect to the temperature evolution, good results 
are obtained with a function that rises rapidly from low temper¬ 
atures to a maximum of just above 30 000 K and then gradually 
drops to about 6 000 K at the end of the flare. Both Jt(t) (right- 
hand scale; normalized to a maximum value of 1; broken line) 
and 0(f) (left-hand scale; solid line) are shown as a function of 
time (the flare duration being normalized to ffl = 1) in the upper 
frame of Fig. 0 

With this choice of temperature and area evolution, and us¬ 
ing the transmission functions of Bessell ( 119901 1. F(t ) (i.e. the 
flare light curve) was calculated for the passbands of the UBVRI 
system. The constant was chosen such that F(t ) obtains a max¬ 
imum of F(t) - 1 in the U band. To mimic a constant back¬ 
ground light, a quantity F^ = 1 was added to all light curves 
before transforming them into magnitudes with an arbitrary zero 
point equal to the magnitude of the first point in the light curve 

3 I chose the symbol © for the temperature to avoid confusion with 
the symbol T that is used in the quantity AT max introduced earlier for 
the time lag of the flickering in different bands. 



Fig. 5. Top: Assumed temperature (solid line, left-hand scale) 
and emitting area (broken line, right-hand scale) evolution of a 
black body used in simulations of the brightness development of 
a flickering flare as a function of time. Bottom: Calculated pro¬ 
files of a simulated flickering flare in the bands of the UBVRI 
system. The inset shows the times of maxima (dots) as a func¬ 
tion of the isophotal wavelength of the respective bands, together 
with a linear least-square fit (broken line). See text for details. 

corresponding to each passband. These light curves, depicting 
the magnitude development of the flare, are shown in the lower 
frame of Fig. [5] 

The simulated flare has several characteristics in common 
with real flickering flares: (i) Bruch (119921) found the mean ratio 
of the gradients of the rise and decline of flickering flares in the 
B band of many light curves (measured between the base and 
the peak of the flare) to be 1.18 + 0.07, meaning that on average 
the rise is slightly more rapid than the decline. In the simulated 
flare a very similar gradient ratio of 1.21 is measured. But note 
that this value depends quite strongly on the choice of ©(f) and 
Jt(t). (ii) The amplitude of the flare is highest in the U band 
and declines monotonically in longer wavelength bands0. This 
agrees with observations. However, it is more difficult to make a 
quantitative assessment here because the amplitude depends on 
the amount of the constant background light, which was (unreal¬ 
istically) assumed to be the same in all passbands in the current 
simulations, (iii) There is a clear dependence of the time of max¬ 
imum on the passband. This is shown in the inset in the lower 
frame of Fig. 0 where the maximum times (dots) are plotted as 
a function of the isophotal wavelength of the photometric band. 
The broken line is a linear least-squares fit to the data. Thus, 
there is a time lag of the flickering observed at different wave¬ 
lengths. The least-squares fit yields dAT^x/d/l = 1.4x 10 -5 time 
units/A. To estimate ffl , the average width cr of a Gauss fit to the 
central peak of the auto correlation function of the B band of 
the light curves of TT Ari used to derive the results quoted in 
Sect. 14.21 was measured: cr = 145 s + 23 s . Taking this as the typ- 
cial flare duration (i.e. ffl = cr), dAr max /dT = 2.1 x 1CT 3 sec/A. 

6 The similarity of the amplitudes in V and R is a consequence of the 
significantly larger width of the R passband. 
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This is of the same order of magnitude as the observed time lag 
forTTAri (Tabled. 

Thus, the simple scenario explored here can explain impor¬ 
tant characteristics of the flickering, in particular the main result 
of the present study, namely the observed time lag of the flick¬ 
ering at different optical wavelengths. The numerical value of 
dATtnax/d/l is of the same order of magnitude as that measured 
in real light curves and depends on the choice of the temperature 
and emission area evolution. Although not investigated here, it 
may be expected that deviations from black-body radiation have 
similar effects. 

The results of this exercise hold for a single flare. But what 
about an entire light curve composed of many flares? If flicker¬ 
ing were to consist of the additive superposition of many inde¬ 
pendent flares, it would be easy to construct such a light curve 
and investigate its properties. However, there are strong indica¬ 
tions in CVs (e.g. Scaringi et al. !2012bt and even more so in X- 
ray binaries (e.g. Uttley et al. ;2005t that flickering is the result 
not of an additive, but of a multiplicative process. The origins of 
disturbances in CV system that gives rise to flickering are largely 
speculative, and the physical mechanisms of how they lead to the 
emission of optical radiation are unknown. Therefore, the prop¬ 
erties of an ensemble of (nonlinearly) superposed flares are not 
straightforward to predict, and I refrain from investigating this 
issue further. 

5.2. Fireballs and the case of SS Cyg 

The heuristic model discussed in the previous subsection can 
obviously not substitute a physical model for flickering flares. 
It was meant to show that parameter combinations for a light 
source exist that can plausibly explain the observed time lag of 
the flickering as a function of wavelength. The evolution of the 
(black body) temperature and the emitting area during the flare 
event was adjusted such that the flicker shape resembles that ob¬ 
served in real light curves. 

A more sophisticated model must take into account the phys¬ 
ical conditions in the flare light source. An interesting attempt 
in this direction has been published by Pearson et al. (l2005l l. 
They modelled the light source as a fireball, that is, a hot, spher¬ 
ically symmetric expanding ball of gas, and calculated the emit¬ 
ting spectrum as a function of the input parameters and of time, 
permitting them to construct light curves at various wavelengths. 
They compared the results with light curves (time resolution: 2 s ) 
in six bands between 3590 A and 7550 A of a particular flare 
extracted from high-speed spectro-photometric observations of 
the dwarf nova SS Cyg. Adjusting the model parameters such 
that the observed light curves are well fitted yields encouraging 
results if the fireball is assumed to be in isothermal expansion. 

To investigate how the results of Pearson et al. (120051 1 fit in 
with the present study, I extracted the light curves of SS Cyg 
(consisting just of a single flare) from their Fig. 14. For conve¬ 
nience, the one corresponding to the band centred on 4545 A is 
shown in the left panel of Fig. [6] In the same way as has been 
done for the other systems investigated in this study, the light 
curves of SS Cyg were cross-correlated, using the one centred 
on 4225 A as reference. The time lag A7 max with respect to the 
reference band is shown as a function of wavelength in the right 
panel of Fig. [6] (filled circles). The wavelength of the reference 
band is marked by a cross. 

For comparison, the average A7 max values for SS Cyg from 
Table [T] are shown in the figure as small dots with error bars 
(mean errors of the mean). The triangle indicates the location of 



Fig. 6. Left frame: A flickering flare of SS Cyg centred on 
4545 A on 1998, luly 8, as observed by Pearson et al. (120051 . 
Right frame: A7 max as a function of wavelength of the compar¬ 
ison band for the observational data (filled circles) and the pre¬ 
ferred fireball model (open circles). The solid and broken lines 
represent a linear least-squares fit to the data points (shortest 
wavelength point excluded). The location of the reference band 
is indicated by a cross. The dots with error bars are the average 
A7 max values of SS Cyg taken from Table Q] with the triangle 
representing the location of the respective reference band. 

the respective reference band. Unfortunately, Table Q] only con¬ 
tains data of SS Cyg taken in the UBV system, defining thus but 
two points in the diagram. Therefore, no strong conclusions can 
be drawn from the comparison with the results extracted from 
Fig. 14 of Pearson et al. (120051 1. However, both data sets appear 
to be roughly compatible with each other. 

Except for the shortest wavelength band centred on 3615 A, 
the data points from Pearson et al. (120051 1 follow a linear rela¬ 
tionship (solid line) with an inclination of (2.3 + 0.1) 10 -3 sec/A. 
This is very similar to the behaviour of V603 Aql, TT Ari, and 
RS Oph (Tableland Fig.©. 

I also determined the times of maximum light of the pre¬ 
ferred model of Pearson et al. (l2005t) by fitting a high-order 
polynomial to the peak of the model curves extracted from their 
Fig. 14. The corresponding A7 raax values are shown as open cir¬ 
cles in Fig. [6] They follow a linear relationship (disregarding 
the data point at 3615 A; broken line) with an inclination of 
(4.86 + 0.01) 10~ 3 sec/A. This is twice the value found for the ob¬ 
servational data, indicating a slight disagreement between model 
and observations. 

Although the fireball model is thus quite attractive to explain 
the time lag of the flickering as a function of wavelength, a vex¬ 
ing question remains: Why does the data point at Pearson et al.’s 
(120051) 3615 A band deviate from the linear relationship between 
A7 nlax and A? Or, more generally, how can the tendency for a de¬ 
viation from such a relationship at wavelengths shorter than the 
Balmer limit (Fig. [3]) be explained? This can be due to optical- 
depth effects in the flaring light source. Pearson et al. (120051) 
found that their fireball is initially dominated by flux coming 
from optically thick regions. This changes in the course of its 
evolution, so that at late times, optically thin emission domi¬ 
nates. Thus, more Balmer continuum emission can escape from 
the light source at later phases, strengthening the flare at ul¬ 
traviolet wavelengths and reversing the AT mllx - A relationship 
observed at longer wavelengths. The balance between optically 
thick and thin radiation will depend on the detailed structure of 
the fireball. Therefore, the late dominance of the Balmer contin¬ 
uum emission may be seen in some, but not all cases, explaining 
the different behaviour of the flickering time lag at very short 
wavelength in individual systems (e.g. V603 Aql and TT Ari vs. 
RS Oph and SS Cyg). 
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5.3. Alternative scenarios 

Scaringi et al. (12013b sought to explain the time lags along lines 
different from fireballs (Sect. or an evolving black-body- 
like light source, as investigated in Sect. 15.II They discarded ex¬ 
planations in the context of models of flickering caused by the 
inward propagation of perturbations within the accretions disk 
[see e.g. Lyubarskii (l!997t >. Pavlidou et al. (12001 b . Kotov et al. 
(120011) and Arevalo & Uttley (120061) for detailed models of this 
kind for X-ray binaries and galactic nuclei, and Yonehara et al. 
(119971 ). Dobrotka et al. (120101 ) and Scaringi (120141) for CVs], 
which assume that the observed optical variations are directly 
related to the inward propagation of matter in the accretion disk 
because then flickering variations in the blue should lag those in 
the red. This is contrary to what is observed. They also rejected 
a scenario of (practically) instantaneous reprocessing of light of 
a variable continuum source close to the centre of the accretion 
disk by some structure farther out. The time lag should then ba¬ 
sically be equal to the light travel time. However, the typical 
dimensions of CVs are such that only the smallest observed time 
lags would not exceed the light travel time from the centre to 
the periphery of the disk. Moreover, the emission at long wave¬ 
lengths should then be dominated by reprocessed light, which 
may not be likely. 

If, however, reprocessing does not occur instantaneously, but 
on the local thermal time scale Scaringi et al. (120131 ) are able 
to qualitatively reconcile the observed time lags with reprocess¬ 
ing sites not too far out, provided that only the surface layer of 
the accretion disk reprocesses photons from the variable light 
source. As a final scenario, Scaringi et al. (120131) mentioned re¬ 
verse shocks in the accretions disk (Krauland et al. 120131 ) that 
may originated in the boundary layer and deposit energy first 
in the hotter inner parts of the disk and then in the cooler outer 
regions. 

None of scenarios mentioned in the previous paragraph is 
backed by a more detailed physical model. Therefore, it is hard 
to say if they lead to time scales compatible with the observed 
lags and are in accordance with other characteristics of the flick¬ 
ering such as, for example, the temporal and spectral develop¬ 
ment of individual flares. But it would be premature to discard 
them before an attempt has been made to predict the implied 
emission properties in some detail and to compare them with ob¬ 
servations. It appears that the fireball model is at present the only 
physical model that has been advanced to a stage where definite 
predictions can be made about at least some flickering proper¬ 
ties and that realistic parameter combinations can be found such 
that these predictions comply with observations. It must be rec¬ 
ognized, however, that the model does not explain what makes 
the fireball explode in the first place. Moreover, it is not known 
to which degree a (linear or nonlinear) superposition of indi¬ 
vidual fireball events can reproduce the statistical properties of 
observed flickering light curves. Thus, it is by no means clear, 
which of the mentioned scenarios, if any, is realized in nature. 

6. Conclusions 

Confirming earlier results of Scaringi et al. (120131 ). it has been 
shown conclusively in the present study that flickering in cata¬ 
clysmic variables does not occur exactly simultaneously in dif¬ 
ferent photometric bands across the optical range. Instead, there 
is a time delay in the sense that individual events of the flickering 
develop slightly later at red than at blue wavelengths. In many of 
the investigated CVs, this effect can only be seen statistically 
in the ensemble of the available data because most of the light 


curves used are not of sufficient quality to resolve the time de¬ 
lay well enough in particular data sets. However, in the systems 
with the best data, V603 Aql and TT Ari, the effect can convinc¬ 
ingly be measured individually. This is also true for RS Oph and 
SS Cyg. To these systems MV Lyr and LU Cam, observed by 
Scaringi et al. (120131 . can be added. 

The time lag of the flickering as a function of wavelengths is 
of the order of a couple of milliseconds per Angstrom. The exact 
value may depend on the particular system, but appears not to 
change much within the same star on time scales of several years 
(i.e. the period spanned by the observational data; this does, of 
course, not preclude time lag changes when systems like TT Ari 
go into a low state). 

The observed effect can be explained if the evolution of the 
radiation characteristics of the light source(s) responsible for the 
flickering is such that individual flares reach their peak emission 
slightly earlier at blue than at red wavelengths. It is not diffi¬ 
cult to construct simple scenarios that are able to reproduce the 
observed time lag not only qualitatively, but also quantitatively, 
and which also agree with other properties that are generally ob¬ 
served in the flickering. While such scenarios lead to the ob¬ 
served effects, they can by no means substitute a realistic phys¬ 
ical mechanism. The fireball model of Pearson et al. (120051) is 
most advanced when it comes to the actual radiation emission 
process that may be seen as flickering. It can reasonably well 
reproduce a flickering flare observed in SS Cyg in various wave¬ 
lengths ranges not only with respect to the wavelength depen¬ 
dent time lag, but also with respect to the spectral and temporal 
development. However, the detailed emission processes implied 
in the context of other scenarios deserve to be investigated as 
well. This hold true especially for models where, contrary to the 
fireball model, a better understanding of the origin of underly¬ 
ing disturbance exists, for instance, the fluctuating disk model. 
It may be interesting to investigate whether the basic fireball ra¬ 
diation physics can be applied to these models as well. 
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